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Ch ■ Abstract 

We investigate the spatiotemporal dynamics of a lattice of coupled chaotic maps whose coupling 

connections are dynamically rewired to random sites with probability p, namely at any instance of 

> 

(•^ ' time, with probability p a regular link is switched to a random one. In a range of weak coupling, 

00 ■ where spatiotemporal chaos exists for regular lattices (i.e. for p = 0), we find that p > yields 



o 



o 



synchronized periodic orbits. Further we observe that this regularity occurs over a window of p 



o 

00 \ values, beyond which the basin of attraction of the synchronized cycle shrinks to zero. Thus we 



have evidence of an optimal range of randomness in coupling connections, where spatiotemporal 



^ . regularity is efficiently obtained. This is in contrast to the commonly observed monotonic increase 

H : 

5r , of synchronization with increasing p, as seen for instance, in the strong coupling regime of the very 

same system. 
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networks, do not to fit in either paradigm 



such as the small-world network 



a. 



INTRODUCTION 

The dynamics of spatially extended systems has been a focus of intense research activity in 
the past two decades. In recent years it has become evident that modelling large interactive 
systems by finite dimensional lattices on one hand, and fully random networks on the other, 
is inadequate, as various networks, ranging from collaborations of scientists to metabolic 

l|. Some alternate scenarios have been suggested. 
Here one starts with a regular structure on a lattice, 
for instance nearest neighbour interactions. Then each regular link from a site is rewired 
randomly with probability p. This model is proposed to mimic real life situations in which 
non-local connections exist along with predominantly local connections. 

There is much evidence that random nonlocal connections, even in a small fraction, sig- 
nificantly affect geometrical properties, like characteristic path length [Sj. However its im- 
plications for dynamical properties is still unclear and even conflicting. While the dynamics 
of coupled oscillators and coupled maps on regular lattices has been extensively investigated 
yl, there have been far fewer studies on the sptiotemporal features of nonlinear elements on 
more general network topologies [5|, 16|] . 

Most existing case studies of coupled networks of dynamical elements indicate that fea- 
tures, such as degree of synchronization, vary monotonically with p. That is, it is observed 
that most dynamical properties interpolate between the limits of regular and random con- 
nections without in any sense being "optimal" or more pronounced at some intermediate 
value of p. 

In this paper however we will provide evidence of a system where there exists a window 
of randomness where one obtains special dynamical features which cannot be anticipated 
from a simple interpolation between the regular and random limits. In particular we will 
show the pronounced enhancement of spatiotemporal order in the system, in am intermediate 
window of rewiring probability p. Our observations are markedly distinct from the commonly 
observed monotonic dependence of synchronization properties on p, as seen for instance, in 
the strong coupling regime of the very same system 6|. 



MODEL 

Specifically we consider a one-dimensional ring of coupled logistic maps. The sites are 
denoted by integers i = 1,. . . ,N, where A^ is the linear size of the lattice. On each site 
is defined a continuous state variable denoted by x„(i), which corresponds to the physical 
variable of interest. The evolution of this lattice, under standard interactions with the 
nearest neighbours on either side, in discrete time n, is given by 

Xn+l{i) = (1 - ^)f{Xn{i)) + -{Xn{i + 1) + Xn{i - 1)} (l) 

The strength of coupling is given by e. The local on-site map is chosen to be the fully chaotic 
logistic map: f{x) = 4x(l — x), as this map has widespread relevance as a prototype of low 
dimensional chaos. 

Now we will consider the above system with its coupling connections rewired randomly in 
varying degrees, and try to determine what dynamical properties are significantly affected 
by the way connections are made between elements. In our study, at every update we will 
connect a site with probability p to randomly chosen sites, and with probability {1 — p) to 
nearest neighbours, as in Eqn. 1. That is, at every instant a fraction p of randomly chosen 
nearest neighbour links are replaced by random links. The case oi p = corresponds to 
the usual nearest neighbour interaction, while p = 1, corresponds to completely random 
coupling. Note that the random connections are dynamic here, as the random links are 
switched around at every update. This is in contrast to most studies which consider static 
(quenched) random connections. 

EMERGENCE OF SYNCRHONIZED CYCLES FROM SPATIOTEMPORAL 
CHAOS 

We will now present numerical evidence that random rewiring has a pronounced effect on 
spatiotemporal order. The numerical results here have been obtained by sampling a large 
set of random initial conditions (~ 10^), and with lattice sizes ranging from 50 to 1000. 

Fig. 1 displays the spatiotemporal state of the network, Xn{i),i = 1, . . .N, with respect 
to coupling strength e, for the case of regular nearest neighbour interactions (i.e. p = 0) 
and for the case of random coupling with probability p = 0.19 and p = 0.6. It is clearly 
seen that the standard nearest neighbour coupling does not yield regularity anywhere in 



the entire coupling range, while randomly rewiring with small probability p = 0.15 creates 
a window in parameter space where synchronized cycles gains stability. Note that different 
periodcities of the synchronized cycles are obatined under different coupling strengths, in the 
window of synchronization. Interestingly however, when p is large this window of complete 
spatiotemporal regularity is lost again (see Fig. Ic). 

In order to quantify this phenomena, we find the fraction of random initial conditions 
that get synchronized (after long transience) This provides a measure of the size of the basin 
of attraction of the synchronized state, denoted by B (see Fig. 2). When B = 1 we obtain 
a global attractor for the synchronized cycles. 

Fig. 3 shows the synchronized basin size B with respect to rewiring probability p, for 
different coupling strengths. It is clearly evident from the figure that for the regular lattice 
limit of p = 0, the synchronized basin size is close to zero. As p > the basin for the 
synchronized state increases rapidly. For instance, for p = 0.01 in Fig. 3, the synchronized 
state is the global attractor of the dynamics, i.e. B = 1. However, as p exceeds an optimal 
window, the basin size decreases rapidly, becoming zero again for high p. Thus there exists 
a range of p values where synchronized cycles gain stability. This observation is quantified 
in Fig. 4 which gives the range of rewiring probability, R, over which B = 1 ior different 
coupling strengths. 

It is evident then that there is a non-monotonic enhancement of synchronization with 
increasing randomness in coupling connections, and synchronization is most enhanced in 
some window of p values. This is in sharp contrast to the monotonic increase in synchro- 
nization with increasing randomness in connectivity, in the strong coupling regime in the 
same system. 

EFFICIENCY OF SYNCHRONIZATION 

We calculate the average time taken for systems with random initial conditions to syn- 
chronize, and denote this as < T >. This quantity, which is a measure of the efficiency of 
synchronization, is displyed in Figs. 5-6. 

It is evident from Fig. 5 that < T > first decreases slightly as p increases and then very 
sharply rises, by orders of magnitude, especially for larger lattices. It is clear then that 
synchronization is most efficient in some range of low p values. So in order to obtain fast 
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FIG. 1: Bifurcation diagram showing the state of the lattice x„(j),i = 1,... A^ {N = 50), over 
n = 1, . . . 25, with respect to couphng strength e, for (a) p = 0, (b) p = 0.19 and (c) p = 0.6. The 
synchronized 2 and 4 cycles are clearly evident gn (b), but not in (a) or (c). 
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FIG. 2: Size of the basin of attraction of the synchronized cycles as a density plot in the space of 
coupling strength e and random rewiring probability p, for a system of size 100. Here transience is 
10000. 

and reliable synchronization it is a good strategy to choose low p, as this ensures a global 
attractor for the synchronized cycles, as well as a fast approach to the synchronized state 
from a random initial state. 

As a consequence of the above, the size of the basin of attraction for the synchronized 
cycles, B, depends on the allowed transience. Fig. 6 shows the increase in the basin B 
as transience is increased. However note that this increase is very slow, and increasing 
transience by 2 orders of magnitude changes the global attractor width in p space by a small 
amount. 

Caveat: it appears possible, from the above observations, that in the limit of T — ;► oo a 
synchronized state may be reached for larger lattices as well. But it should be underscored 
that transience is always finite, and it is indeed very pertinent to know what generically 
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FIG. 3: Size of the basin of attraction of the synchronized state vs rewiring probabihty p for a 
system of size 100 for coupHng strengths (a) e = 0.16, (b) e = 0.19 (c) e = 0.20. Here transience is 
5 X 10^. 

happens in the limit of finite, but large (~ 0{N'^)), transience, as is the case for the studies 
in this paper. 

CONCLUSIONS 



We have investigated the spatiotemporal dynamics of a network of coupled chaotic logis- 
tic maps, with varying degrees of randomness in coupling connections. We find a window 
in coupling parameter space, in the weak coupling regime, where random rewiring induces 
spatiotemporal order. Interestingly the basin of attraction for the synchronized cycles varies 
non-monotonically with rewiring probability p. As p is increased, the basin of attraction of 
the synchronized state rapidly increases. At an intermediate (small) value of p the synchro- 
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FIG. 4: Range of rewiring probability R, over which synchronization is achieved vs. coupling 
strength e. Here system size is 100 and transience is 10000. 

nized state becomes the global attractor of the system with (almost) all initial conditions 
being attracted to the synchronized state. However, interestingly, as p is further increased, 
the basin of attraction of the synchronized state shrinks to zero again. Thus we have strong 
evidence of the pronounced enhancement of spatiotemporal order in an intermediate window 
of p. This is in marked contrast to the monotonic increase in synchronization with increasing 
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FIG. 5: Average time required to reach synchronization < T > vs rewiring probabilty p, for 
couphng strength e = 0.19 for system sizes N = 25,50,100. The synchronized state is a gfobal 
attractor for the system for the values of p displayed here. 

randomness in connectivity, observed in the strong coupling regime of the same system. 
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FIG. 6: Size of the basin of attraction of the synchronized state B vs rewiring probabihty p, for 
couphng strength e for a system of size 100 for maximum allowed transience equal to (a) 5 x 10^ 
(b) 5 X 10^ (c) 5 X 10^. 
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